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Abstract 



The production of prompt J/V', i^C^S), Xci and Xc2 is studied using a 
32.4 fb~^ data sample collected with the Belle detector at the T(45') and 
60 MeV below the resonance. The yield of prompt J/ijj mesons in the T(4S') 
sample is compatible with that of continuum production; wc set an upper 
limit B{T{4S) J/i^X) < 1.9 x 10"^ at the 95% confidence level, and find 
a{e^e~ — > J/ip X) = 1.47 ± 0.10 ± 0.13 pb. The cross-sections for prompt 
^{2S) and direct J/ip are measured. The J/ip momentum spectrum, produc- 
tion angle distribution and polarization are studied. 

PACS numbers: 13.65.+i, 13.25.Gv, 14.40.Gx 
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Typeset using REVT^iX 



The production of prompt charmonia is poorly understood, and provides an interesting 
environment to study the interplay between perturbative QCD and non-perturbative effects. 
A recently developed effective field theory, called non-relativistic QCD (NRQCD) [Q, pro- 
vides a consistent calculational framework for direct heavy quarkonium production. Further 
experimental information is needed to establish the applicability of NRQCD to charmonium 
production 0. Studies of prompt charmonia in e+e~ collisions at the T{4S) resonance con- 
stitute a test of NRQCD and can provide estimates of some of its non-perturbative matrix 
elements. 

In this Letter, we report results of a measurement of prompt charmonium production 
in data recorded at the T(45') and in the continuum 60 MeV below the resonance, corre- 
sponding to integrated luminosities of 29.4 fb~^ and 3.0 fb~^ respectively [^f. The data was 
collected with the Belle detector at the KEKB asymmetric energy (3.5 GeV x 8 GeV) e^e~ 
storage ring [Q. 

The Belle detector is a large-solid- angle spectrometer equipped with a 1.5 T super- 
conducting solenoid magnet. Charged tracks are reconstructed in a 50 layer Central Drift 
Chamber (CDC) and in three concentric layers of double sided silicon strip detectors (SVD). 
Photons and electrons are identified using a CsI(Tl) Electromagnetic Calorimeter (ECL) lo- 
cated inside the magnet coil. Muons and mesons are detected using resistive plate 
chambers embedded in the iron magnetic flux return (KLM). Charged particles are iden- 
tified using specific ionization measurements in the CDC, pulse heights from the Aerogel 
Cerenkov Counters (ACC) and timing information from the Time of Flight Counters (TOF). 
A detailed description of the Belle detector can be found elsewhere [|]. 

Hadronic events are separated from QED, rr, two-photon and beam-gas interaction 
backgrounds by requiring the presence of at least three charged tracks (A^^ch > 3), an event 
vertex with radial (r0) and z coordinates within 1.5 and 3.5 cm of the origin, respectively, a 
total reconstructed center-of-mass (CM) energy greater than 0.2^/s {^/s is the CM collision 
energy), a z component of the net reconstructed CM momentum less than 0.5y^/c, a total 
ECL energy between O.l^i and 0.8a/s with at least two energy clusters associated, and i?2, 
the ratio of second and zeroth Fox- Wolfram moments, less than 0.8. 

Candidate J/ip mesons are reconstructed using the leptonic decays J/ip ^ jji^ and 
e^e~ . For J/ip f^^f^~, both charged tracks must be identified as muons in the KLM 
using information on hit positions and penetration depth. For J/ip ^ e~^e~, oppositely 
charged track pairs must be identified as electrons based on a combination of CDC dE/dx 
information, ACC response, and the position, shape and energy of the associated ECL 
shower. To correct for final state radiation and bremsstrahlung, photons within 50 mrad 
of the are included in the e^e~ invariant mass calculation. The two lepton candidate 
tracks are required to have a common vertex, with a distance in the r0 plane to the average 
interaction point dr,j, < 500 /xm. The signal region is defined by the mass window —93 < 
Mi+i- — Mjj^ < 33 MeV/c^, common for both decay channels. 

We reconstruct ip{2S) J/ipn^TT" decays by combining J/ip candidates with tt+tt^ pairs 
and requiring a common tt+tt^ vertex. Candidate Xci and Xc2 mesons are reconstructed via 
their radiative decays to the J/ip: we combine J/ip candidates with photons detected in the 
ECL, that are not associated with an identified 7r°. 

The largest source of secondary charmonia in the T(4S') sample, due to B meson decays, 
is eliminated by requiring the charmonium CM momentum p* to be above the B decay 
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kinematic limit. A common requirement p* > 2.0 GeV/c is used for all analyses {J/ip, V'(2S'), 
Xci and Xc2)'- this value is robust against the effects of momentum measurement errors and 
motion of the S-meson in the CM. This requirement is not applied to off-resonance data. 

The background due to initial state radiation with a hard photon ("radiative return to 
J/ip {iIj{2S))") P and higher order QED processes e~^e^ J/ip-y*, J/ipl^l^ [0] is also large. 
This background contributes mainly to Nch = 3 and Nch = 4 events, and a dedicated study 
was performed for each of these samples. The J/iJj mass peak for Nch = 3 has a signal to 
background ratio S/ B 0.1, which is too poor to be included in the final signal sample. For 
A^ch = 4 the S/ B ratio is acceptable, but 85% of the J /ip mesons are due to QED processes, 
principally e^e~ — > tp{2S)'~f —>■ J/tpTc^iT^ (7) with or without a visible photon. Based on 
this study we require Nch > 4 to suppress the QED backgrounds, and account for the loss of 
low multiplicity signal events in the efficiency calculation. The remaining QED background 
in Nch > 4 events is suppressed by rejecting events with a 'ipi^S) J/ip7i~^n~ candidate 
accompanied by either a photon with CM energy E* > 3.5 GeV or electron tracks from 
photon conversion. 

Mass distributions for prompt J ftp candidates in the T(4S') and continuum data samples 
are shown in Fig. |l|. Clear signals are observed both in the off-resonance (226 ± 23 events) 
and the T(4S') data (1459 ± 57 events), for both decay channels. The average detection 
efficiency for J/^p e+e" (fi+f^-) is 38% (48%). 

The ip{2S) signal in T(45') data with p* > 2 GeV/c is shown in Fig. ^(a). Combining 
J/ip ^ e~^e~ and clear signal of 143 ±19 events is seen in the mass difference 

M/+;-^+7r- — Mi+i- distribution. In the off-resonance data, no statistically significant ip{2S) 
signal is seen (10 ± 5 events); this is not inconsistent with the observed rate at the T(4S') 
resonance. Searching for Xci and Xc2 in the T(4S') data, we form the mass difference Mi+i-^ — 
M1+1-] no significant signals are seen (Fig. 0(b)). The average detection efficiency for il){2S) 
and Xci, Xc2 is 20% and 21% respectively. 

We first compare J/ip production for p* > 2 GeV/c in the T(45') and continuum data 
samples by extracting the number of signal events from invariant mass distributions for 
300 MeV/c wide p* bins. The shape and normalization of the resulting p* distribution for 
the off-resonance continuum (not shown) agree very well with those for the T{AS) data, 
indicating that all the J/ip production at the T(45') can be explained by continuum pro- 
duction. The branching fraction for the decay T(4S') J/ipX can be estimated from the 
difference between the normalized yields (1459 ± 57 and 1496 ± 145 from the T(4S') and off- 
resonance data respectively) which is found to be — 37±156. We therefore set an upper limit 
B{T{AS) J/^X) < 1.9 X 10-^ (at 95% CL.) using the Feldman-Cousins method 0. This 
limit is more stringent than the result of Ref . . There are no published predictions for 
B{T{AS) J/ipX), but scaling the NRQCD estimate B{T{IS) J/^X) ~ 4.5 x 10"^ |ll| 
by the ratio of widths rx(is)/rT(4S) gives ~ 1.7 x 10~^; this is far below our current experi- 
mental sensitivity. 

Hereafter we assume that all prompt J/ip's produced at ^/s ~ 10.6 GeV are due to con- 
tinuum production. We combine off- resonance data with the T(4S') data for p* > 2 GeV/c 
to extract charmonium production cross-sections and angular distributions. Acceptance 
corrections are determined using a dedicated Monte Carlo (MC) simulation of the process 
e~^e~ — > J/ip {■iIj{2S)) qq, where the composition of qq flavors is set to that observed at 
^/s ~ 10.6 GeV. To minimize model dependence of the efficiency determination, the data is 
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corrected using a two-dimensional acceptance weight matrix as a function of p* and cos^*, 
where 9* is the J/ip production angle in the CM system. The momentum range p* > 2 GeV/c 
is divided into three bins, and cos 9* into five bins; the number of signal events in each two- 
dimensional bin is obtained by fitting the invariant mass distribution. A similar procedure is 
applied to determine the ip{2S) acceptance for p* > 2 GeV/c, while acceptance corrections 
for Xci,c2 were estimated directly from the MC without binning. 

The resulting cross-sections are summarized in Table |I[ For the J/ip sample {p* > 
2 GeV/c), the feed-down from prompt ^{2S) is measured to be 0.33 ± 0.04to;o6 pb, and this 
contribution is subtracted to give the direct J/ if) cross-section in the table. The p* <2 GeV/c 
off-resonance data is used to extend the J /ip cross-section measurement to the full p* range. 
As we observe no significant signal in the continuum for ip{2S)^Xci and Xc2, we estimate 
cross-sections for p* > 2 GeV/c only. 

For the systematic errors of the measured cross-sections the following contributions are 
considered (where appropriate): uncertainty in the efficiency determination (±5%), effects of 
the multiplicity cut (+4%), lepton identification (±4%), tracking (±4%), luminosity mea- 
surement (±1.3%), possible feed-down from the Xci,c2 (—10%), contamination from QED 
processes (—4%), ip{2S) feed-down subtraction (±2%) and errors of the respective branch- 
ing fractions. A contribution from the unobserved process T(4S') —>■ J/ipX is not included. 
If the true branching fraction for this decay were just below the determined upper limit, 
the measured J/ip cross-sections would be overestimated by ~ 0.18 pb. Other sources of 
systematic error, such as contamination from radiative return to the T(1S', 25", 35") and from 
77 ~^ Xc2, were found to be negligible. As an additional check, J/ip cross-sections were 
determined separately for the e~^e~ and decays, and found to be in good agreement: 

their ratio is 1.00 ± 0.07 ± 0.04. 

Our results on the J/ip total and partial cross-section are smaller than those of 
BaBaR ||10|. The partial cross-sections for direct J/ip and prompt ip{2S) production rep- 
resent the first measurements in e~^e~ collisions. 

NRQCD predictions for a{e~^e~ J/ipdhect^) cover the range 0.8 — 1.7 pb ||12[ |13 



ri| . (We exclude the extreme value of Ref. [0.) The most recent analysis |11| obtains 



a total cross-section of 0.8 — 1.1 pb, attributing 0.3 pb and 0.5 — 0.8 pb to color-singlet 
and color-octet mechanisms respectively. Our result can be used to further constrain the 



combination of color-octet matrix elements needed in the calculation [jT2[ |]TT|. The ratio 
a{'ilj{2S) X)/a{J/ipdirectX) is measured to be 0.93 ± O.lltfl;}^; this is related in NRQCD 
to a ratio of linear combinations of non-perturbative matrix elements for the respective 
mesons 
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Momentum distributions are shown in Fig. ^ The distribution of feed-down from ip{2S) 
to J/ip is also shown in Fig. |5(a): it does not significantly modify the overall shape of 
the J/tp momentum distribution. The J/ip distribution vanishes ~ 300 MeV/c below the 
kinematical limit (4.84 GeV/c) while the ip{2S) momentum distribution extends till the 
end-point (4.65 GeV/c). The J/tp distribution is softer than the NRQCD prediction for 
color-singlet J/ip gg ||T3[ [|T^, and agrees quahtatively with the predicted shape of the color- 



singlet J/ip cc component |jT5| |]T3|. Some NRQCD calculations also predict a dramatic rise in 
the cross-section at the end-point due to color-octet e~^e~ J /ip g [0] |]I5| . The detection 



efficiency for this process was studied using a dedicated generator embedded in Pythia [17 



Assuming a 1 pb cross-section for this process, we expect > 300 events in the last two bins 
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of Fig. 15(a). No such signal is observed. 

The distributions of the prompt J /ip CM production angle 9*, and the helicity angle 9h 
(the angle between the positive lepton daughter momentum vector in the J /if) rest frame, and 
the J /if) momentum vector in the CM system), have also been studied in different momentum 
intervals. We correct the distributions for detection efficiency, but ignore possible effects 
due to feed-down. The distributions are fitted with the parameterizations 1 + A cos^ 9* and 
1 + a cos^ 9h'i the results are summarized in Table with selected fits also shown in Fig. ^. 
No statistically significant p* dependence is seen for either A nor a parameters. 

A large positive A at all momenta for direct J/ip production is expected only for the 
color-singlet J /if) cc mechanism [jT5[. However its contribution to the cross-section is thought 
to be small (~ 10%) [|l^ [l^. The leading color octet process J/'ipg is expected to yield 
A ^ +1 at the endpoint |T6[ although as noted above we do not observe these events in 



the p* distribution. A significant longitudinal polarization of direct J/ip mesons [a < —0.4) 
is expected for the color-singlet process J/ifj gg alone [|T^]. 

In summary, we have observed production of prompt J /ip and ip{2S) at energies near the 
T(4S') mass. We set an upper limit for J /ip production from T(4S') (valid for p* > 2 GeV/c), 
B{T{AS) J/^X) < 1.9 X 10"^ at the 95% C.L., and find the total cross-section for contin- 
uum prompt J/^/) production to be (T(e+e~ — > J/ipX) = 1.47ib0.10(stat.)±0.13(syst.) pb. In 
the momentum range p* > 2 GeV/c, we measure a{e^e^ J/ipdhect X) = 0.72±0.0Sll'^[\j pb 
and(T(e+e- ^ ij{2S) X) = 0.67±0.09i[]:?? pb, and their ratio (t(^/'(25) X)/a( J/^/^direct ^) = 
0.93 ±0.171 The angular distribution of prompt J/ip mesons follows 1 + Acos"^ 9* with 
A = 0.9 ± 0.2, and the helicity angle follows 1 + acos^ 9h with a = —0.4 ± 0.1 indicating 
partial longitudinal polarization. We do not observe a statistically significant variation of A 
or a with momentum. 

We wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator. We acknowledge support from the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan and the Japan Society for the Promotion of Science; 
the Australian Research Council and the Australian Department of Industry, Science and 
Resources; the Department of Science and Technology of India; the BK21 program of the 
Ministry of Education of Korea and the CHEP SRC program of the Korea Science and 
Engineering Foundation; the Polish State Committee for Scientific Research under contract 
No.2P03B 17017; the Ministry of Science and Technology of Russian Federation; the National 
Science Council and the Ministry of Education of Taiwan; and the U.S. Department of 
Energy. 
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TABLES 



TABLE I. The measured charmonium cross-sections. The upper hmits for Xci and Xc2 are at 



the 90% C.L. 



a 




0<p* < p*^^^ 


2.0 <p*< p*^^ 


[Pb] 




[GeV/c] 


[GeV/c] 


a{e'^e~ 




1.47 ±0.10 ±0.13 


1.05 ± 0.04 ± 0.09 




J/'fpdirectX) 




0.72 ± 0.08l|]:J^ 


a{e^e~ 


^ i^{2S)X) 




0.67 ± omtoil 




^ XciX) 




< 0.35 


o"(e+e~ 


^ Xc2X) 




< 0.66 



TABLE IL Results of the fits to angular distributions. 



P}/^ [GeV/c] 


A 


xVd.o.f. 


a 


xVd.o.f. 


2.0 - 2.6 


n 0+U.5 


1.5/4 


-0.4 ±0.2 


7.8/4 


2.6 - 3.4 


1 1+°-^ 

^•^-0.3 


5.0/4 


-0.4 ±0.1 


1.3/4 


3.4-4.9 


1 1+0.4 
-'■•-'■-0.3 


4.5/4 


-0.2 ±0.2 


7.1/4 


2.0-3.4 


0.7 ±0.3 


1.2/4 


-0.5 ±0.1 


4.6/4 


2.0 - 4.9 


0.9 ±0.2 


3.0/4 


-0.4 ±0.1 


13.4/4 
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FIGURES 



i 1 1 1 i 1 1 
■ Nj,^ = 651 ± 41 (a) 
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. Nj,„ = 808±39 
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■ Nj,^=102±16 (c) . 

1 . . . 


I'll 

■ ^/„=124±16 

1 


(d) . 



2.60 3.00 2.60 3.00 3.40 



M^^^_ (GeV/c") M^^^_ (GeV/c") 

FIG. 1. Mass distributions for the J/^ ^ e^e~ (a,c) and f^^fJ-^ O^A) candidates; (a) and (b) 
are for T(45) data with p* > 2 GeV/c, (c) and (d) are for off-resonance data. Nj/^ is the number 
of J/ip mesons determined by a fit to the dilepton mass distribution, where a Crystall Ball line 
shape function Q is used for the signal and the background is described by a polynomial. 
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Nx„ = 31 ± 24 




Nx,2 = 23 ± 24 J 




■~J 1 


■ , , , , 1 


, , , , 1 , , , , " 



0.60 



FIG. 2. Mass difference distributions for charmonium candidates in the T(4S') data: (a) 
Mi+i-^+^~ - Mi+i- for '>i:{2S) candidates with p* > 2 GeV/c. (b) Mi+i-^ - Mi+i- for Xci,c2 
candidates with p* > 2 GeV/c. The curves represent fit results with a Crystall Ball line shape for 
the signal and a polynomial for the background. 
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FIG. 3. CM momentum distributions of prompt charmonia, corrected for efficiency : (a) J /ip 
(filled points) and J/ip mesons from ip{2S) — > J/ip X (open points); (b) iIj{2S). 
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FIG. 4. The cos^* (a) and cos 9h (b) distributions for low and high p* . The curves represent 
fit results, described in the text and summarized in Table ||. 



11 



